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Abstract 

Traditional quantum error-correcting codes are designed for the depolarizing channel modeled by generalized 
Pauli errors occurring with equal probability. Amplitude damping channels model, in general, the decay process 
of a multilevel atom or energy dissipation of a bosonic system at zero temperature. We discuss quantum error- 
correcting codes adapted to amplitude damping channels for higher dimensional systems (qudits). For multi-level 
atoms, we consider a natural kind of decay process, and for bosonic systems, we consider the qudit amplitude damping 
channel obtained by truncating the Fock basis of the bosonic modes to a certain maximum occupation number. We 
construct families of single-error-correcting quantum codes that can be used for both cases. Our codes have larger 
code dimensions than the previously known single-error-correcting codes of the same lengths. Additionally, we present 
families of multi-error correcting codes for these two channels, as well as generalizations of our construction technique 
to error-correcting codes for the qutrit V and A channels. 

Index Terms 

amplitude damping channel, quantum codes, qudit 


I. Introduction 

For a q-level quantum system with Hilbert space C^, called a qudit, the most general physical operations (or 
quantum channels) allowed by quantum mechanics are completely positive, trace preserving linear maps which can 
be represented in the following Kraus decomposition form M(p) = ^kpE^, where the matrices Ek are called 
Kraus operators of the quantum channel M satisfying the trace-preserving condition £'fc£'fc = 1- 

In designing error-correcting codes for protecting messages carried by n qudits sent through a channel J\f, it is 
usually assumed that the errors to be corrected are completely random, with no knowledge other than that they 
affect different qudits independently [16, 30]. The corresponding channel J\f is the depolarizing channel which can 
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be modeled by uniformly distributed error-operators given by generalized Pauli operators [17, 25, 31] (Xg)°'(Zg)^, 
for a,b € {0,1,..., g — 1}, where Xq\s) = |s -f 1 mod q), Zg\s) = w®|s), and w = exp(27ri/g). When it is clear 
from the context, we may just write X and Z, dropping the index q. 

However, if further information about the error process is available, more efficient codes can be designed. Indeed, 
in many physical systems, the noise is likely to be unbalanced between amplitude (W-type) errors and phase (Z- 
type) errors. Recently a lot of attention has been put into designing codes for this situation and into studying their 
fault tolerance properties [1, 12, 15, 20, 34, 40]. All these results use error models described by Kraus operators 
that are generalized Pauli operators, but for those error models, the W-type errors (i.e., non-diagonal Pauli matrices) 
happen with probability which might be different from the probability p^ that Z-type errors (i.e., diagonal Pauli 
matrices) happen. The quantum channels described by this kind of noise are called asymmetric channels. 

A closer look at the real physical process of amplitude damping noise shows that one needs to go even further, 
beyond Kraus operators of Pauli type. To be more precise, for q = 2, the qubit amplitude damping (AD) channel 
is given by the Kraus operators [9] 



( 1 ) 


Since the error model of the qubit AD channel is not described by Pauli-type Kraus operators, the known techniques 
dealing with Pauli errors result in codes with non-optimal parameters. Several new techniques for the construction 
of codes adapted to this type of noise with non-Pauli Kraus operators, and the qubit AD channel in particular, have 
been developed [9, 15, 27, 28, 36]. Systematic methods for constructing high performance single-error-correcting 
codes [27, 36] and multi-error-correcting codes [11] have been found. 

In this paper, we discuss constructions of quantum codes for AD channels of general qudit systems. Unlike 
the qubit case, where the AD channel is unique, for qudit systems there are different AD channels associated 
with different physical systems. We will focus on two different models: multi-level atoms with a natural kind of 
decay process, and bosonic systems obtained by truncating the Fock basis of the bosonic modes to the maximum 
occupation number q — 1 for a single bosonic mode. 


II. The Amplitude Damping Channel 


For two-level atoms, the decay process at zero temperature is described by the Kraus operators Aq , Ai as given 
in Eq. (1). For multi-level atoms, there are different kinds of decay processes at zero temperature. One natural 
decay process is the cascade structure 5, where the decay process is governed by the master equation [32, 35] 



( 2 ) 


Here {|*)}f^Q is a basis of the Hilbert space C^, cr^_^ i = I* ~ l)(*l ^nd cr^_i j = \i){i — 1|. 


The solution to this master equation gives the Kraus expression 


:=.{p) — AopAl + ^ AijpAlj, 


0<i<j<q—l 


(3) 
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where Aij = y/^\i){j\ with positive coefficients jij, and Aq is a diagonal matrix with its elements determined 
by AlqAq + J2o<i<j<q-i ~ Furthermore, when the decay time t is small, jij is of order for any 

j = i + i, £ > 0. As a consequence, Aq is of order r, and Aij is of order for any j = i + i, £ > 0. This is 
intuitively sound as for the cascade structure, the first order transition always happens from \i + 1) to |i). 

As an example, for three-level atoms, i.e., g = 3, we have 

^01 = Vt^|o)(i|, 

^02 = \/7o2|0)(2|, Ai2 = \/7i2|1)(2|, 

^0 = |0)(0| + \/l — 7oi|l)(l| + \/r^^7o7^^T^|2)(2|, 

where 


7 01 = 2k2T + O(r^), 

702 = 2kik2T^ + O(r^), 
712 = 2 kiT + O(r^), 


for fci 7 ^ /c 2 - The values of jij are slightly different for ki = k 2 , but the order of jij in terms of r remains the 
same. 

The channel A describing energy dissipation of a bosonic system at zero temperature is discussed in [9]. The 
Kraus operators are given by 

9-1 

Ak = J2 

r—k 

where q — 1 is the maximum occupation number of a single bosonic mode, and fc = 0,l,...,g—1. The parameter 
7 is of first order in terms of the decay time r, i.e., 7 = cr + O(t^). As a consequence, the non-identity part of 
Aq is of order r, and Ak is of order for 1 < fc < d — 1. 

For instance, for the qubit case, i.e., g = 2, we have the qubit amplitude channel given by Eq. (1). For g = 3, 
we have 




— k^/k 


ryk\j. _ k){r\, 


(4) 


Ao = |0)(0| + vT^|1)(1| + (1-7)|2)(2|, 

Ai = V7|0)(l| + v/ 27(1-7)|1)(2|, 
and A 2 = 7 | 0 )( 2 |. 

Note that for g = 3, the non-diagonal Kraus operators of the channel A for bosonic systems are linear 
combinations of the Kraus operators of the channel S. Hence codes correcting errors of the channel S are also 
codes for the channel A. 


III. Error Correction Criteria 

A quantum error-correcting code Q is a subspace of (C®)®”, the space of n qudits. For a iC-dimensional code 
space spanned by the orthonormal basis \tjji), i = 1, • ■ •, A', and a set of errors £, there is a physical operation 
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correcting all elements Ek G 8 if the error correction conditions [2, 26] are satisfied; 

yi,j,k,l: {i{ji\ElEi\tj}j) = Sijaki, (5) 

where aki depends only on k and 1. A code is said to be pure with respect to some set of errors £ if au = 0 for 

k A if-dimensional code with length n is denoted by {{n,K)). 

For the AD channels S and A, if the decay time t is small, we would like to correct the leading order errors that 
occur during amplitude damping [3, 4]. Similar as for the qubit case [16, Section 8.7], we will show below that in 
order to improve the fidelity of the transmission through the AD channel S or ^ from 1 — 0 (t) to 1 — 
i.e., to correct t errors, it is sufficient to be able to detect t errors of type Aq and to correct up to t errors of type 
Aij with j > j of total order for the channel S (or to correct up to t errors of type Ai with * > 0 for the 

channel A). We will then say that such a code corrects t amplitude damping errors since it improves the fidelity, 

just as much as a true f-error-correcting code would for the same channel. This is a direct consequence of the 
following sufficient condition for approximate error correction. 

Theorem 1: Assume we are given a quantum channel with Kraus operators Ek that have a series expansion in 
terms of ^/T for a parameter r. A quantum code Q with orthonormal basis {|ci) ■. i = 1,..., K} corrects errors up 
to order 0(r*) if the following conditions are fulfilled for all basis states |ci), \cj) and all error operators Ek,Ei: 

{ci\ElEi\cj) = SijVki + ( 6 ) 

Proof: Assume we are given a quantum channel with Kraus operators Ekir) that depend on some small 
parameter r. We expand the operators in terms of y'r as 

Ekir) = EkmT"^^^- (7) 

m>0 

This leads to the following description of the channel: 

k k m>0 /i.>0 

Given a quantum code Q with orthonormal basis {jci): z = 1,..., K}, the conditions for perfect error correction 
are 

{c^\El{T)Ei{T)\cj) = Sijaki. (9) 

Using the expansion of the Kraus operators in terms of we get 

E = kjakt. (10) 

We are looking for sufficient conditions such that the residual error with respect to p is of order for 


some f > 0. 
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We write each Kraus operator Ek = Bk + Ck + Dk as sum of three terms, where 


Bk=Yl EkmT”^'\ 

m—0 
2t 

Ck= Y. 

Dk= Y 

m>2t 

Then the original channel Af can be written as 

Af{p) = ^ ^ Bk + Ck + Dk)p{Bk + Ck + Dk)^ 

k 

= Y BkpBl 

k 

+ Y^ BkpCl + CkpBl 

k 

+ Y^ BkpD\ + DkpBl 

k 

+ Y^'^k^^k)p{Ck+Dk)^. 

k 

Condition (6) implies that 


( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 


(19) 

( 20 ) 


^^i\BkBl\Cj) ^ijAkl-, 

{ci\BlCi\cj) = SijPki + 0(r*+^). 

In particular, the error operators Bk can be perfectly corrected. We first define the projection operator onto one of 
the spaces' spanned by {i?fc|ci): i = 1 ,... ,K}: 

PB,=YBkh){c.\Bl (21) 

i 

and the partial isometry that maps Bk\ci) to \ci): 

t/B. =^|Q)(c.|i?i. (22) 

i 

We compute 


Ub,Pb, = (^Bk\c,){c,\B 

= \c){ci\B^i^Bk\cj){cjlB^, 


(23) 

(24) 




= XkkY\c){c^\Bl (25) 

i 

The last step follows from the fact that the error operators Bk can be perfectly corrected, which also determines 

^We may, without loss of generality, use lineai' combinations of the original error operators Bk such that these spaces become mutually 
orthogonal. 
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the constant X^k- Then the partial correction operator TZp^a is given by 


For a general state 


^part(p) — ^ ' \Xkk\ 'y ^ P ^k\Cj){Cj\. 

k i,j 


(26) 


in the quantum code Q, the term (15) of M{pq) reads 

J^b^Pq) = 

k i,j 

Since the error operators Bk can be perfectly corrected (implied by Eq. (19)), it can be shown that applying the 
partial recovery operator to J^B^iPo) yields a state XpQ that is proportional to the original state pq. Hence after 
partial recovery we have 


'^part(A/'(pQ)) = XpQ +S{pq), (29) 

where the map S is given by the application of the partial recovery operator to the terms given in (16), (17) and 
(18). The summands (17) and (18) are all of order 0(r*+^), so we can ignore them, but (16) contains terms of 
order for t < I <2t. Applying the partial recovery operator to (16) and using (20) results in the state 


^part i? BkPQCl + CkPqBlj 

= ^^|Afefe| \ci){ci\Bl. I BipqCj + CipqBjj Bk\cj){cj\ 
k i,j \ I / 

— \ ^kk\ i^j' \ ^j ^k\Cj){Cj \ + jc^) (q |1 Cjf ) {Cj/ \ Sj Bk \ Cj ) {Cj \ ^ 


i',j' 




k,l i',j' 


= PQ ^ \^kk'^ {XklP*kl + ^klPkl) + 


k,l 


(30) 


which is, up to order 0 (t*“'"^), proportional to the original state. 

■ 

Note that in the proof of Theorem 1, we have split the error-operators accordingly based on their expansion (7) 
in terms of see (11)-(13). Clearly, the high order parts Dk can be completely ignored. Only the errors Bk of 
approximately half the final order have to be corrected (19), while the errors Ck have to obey some kind of error 
detection criterion (20). 
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IV. Stabilizer and Asymmetric Quantum Codes 

Before presenting our construction of quantum codes tailored to amplitude damping channels, we investigate the 
performance of traditional quantum error-correcting codes on these channels. 

Stabilizer codes are a large class of quantum codes which contain many good quantum codes [16, 30]. A stabilizer 
code Q with n qudits encoding k qudits has distance d if all errors of weight at most d—1 can be detected or have 
no effect on Q, and we denote the parameters of Q by |n, k, djg. Obviously a stabilizer code of distance 2t + 1 
corrects t AD error as it corrects t arbitrary errors. 

Calderbank-Shor-Steane (CSS) codes [6, 37] are a subclass of the stabilizer codes. It has been shown that CSS 
codes can be used to construct codes for the binary AD channel [16, Section 8.7]. The construction is based on 
so-called asymmetric quantum codes, which have a direct generalization to the qudit case [13]. The following 
theorem shows that those asymmetric CSS codes can also be used to obtain error correcting codes for qudit AD 
channels. 

Theorem 2: An |n, kjq CSS code Q with pure -distance 2f-|- 1 and pure Z-distance t corrects t AD error, i.e., 
errors up to order 0 (t‘). 

Proof: The generalized Pauli operators form a basis for all operators on a single qudit. Hence we can 

expand the error operators Ai in terms of tensor products of the generalized Pauli operators. The diagonal error 
operator Aq of AD channels can be expanded in terms of the error operators Z*, with the expansion coefficients 
of the operators Z', / > 0 being of first order in r. The diagonal of the other error operators or Ai is zero. 
They can be expanded in terms of operators Z^X^', fc ^ 0, with the expansion coefficients being of order ^/t. 

Note that for CSS codes, X and Z errors can be corrected independently. The error operators Bk defined in (11) 
of the proof of Theorem 1 are of order at most f/2 in r, and hence they contain no more than t X-errors and no 
more than t/2 Z-errors. As the code Q has X-distance 2f -|- 1 and Z-distance t + 1, the error operators Bk can be 
corrected. Similarly, for the error operators Ck defined in (12), there are no more than 2t X-errors and no more 
than t Z-errors, which can be detected using Q. Hence, using Theorem 1, it follows that Q corrects all errors up 
to order 0 (t*). ■ 


V. Classical Asymmetric Codes 

In this section, we construct quantum codes correcting a single AD error using classical asymmetric codes. Codes 
for the qubit case have been presented in [27, 36]. Those codes are self-complementary, i.e., the basis states are of 
the form Iftf) = '^(1'^) + 1^))’ where u is an n-bit string, u = 1 0 u, and 1 is the all-one string. 

For the non-binary case with g > 2, we consider a similar construction. Define X = X®”, then the basis states 
are chosen as ^ 

IV'u) = ^^^'|u)- 

For instance, for g = 3 and n = 3, we get \tj}o) = -^(lOOO) + |111) + |222)). 


( 31 ) 


The quantum code Q is then spanned by {|'0u)}, where u S C is some length-n string over the alphabet 
{0,1,..., q— 1} (C is a classical code of length n). The advantage of this construction is that the code automatically 
satisfies the error-detection condition for a single Z* error (I = 1,2,... ,q — 1), as the code is stabilized by X. 
Now consider a classical code with codewords C = {u + al: u € C, a = 0,..., q — 1} and the corresponding 
quantum code spanned by {{ipu ): u G C}. The problem of correcting a single error for the qudit AD channels can 
then be reduced to finding certain classical codes. 

The relevant classical channel is the classical asymmetric channel [23]. Let the alphabet be with the ordering 
0<l<2<---<g—l.A channel is called asymmetric if any transmitted symbol a is received as 6 < a. The 
mostly studied asymmetric channel, dating back to Varshamov [38], can be described by the following asymmetric 
distance A(x, y). 

Definition 3 (see [21 ]): Let B = {0,1,..., g — 1} C Z. For x, y G i?", we define 

1) w{x) := Er=i2^*- 

2) iV(x, y) := max{yi - 0}. 

3) A(x,y) := max{A(x,y), A(y,x)}. 

If X is sent and y is received, we say that w(x — y) errors have occurred (note that Xi > yi and hence each 
summand in w(x — y) is nonnegative). A code correcting f-errors is called a t-code. 

Theorem 4 (see [21 ]): A code C C i?" corrects r errors of the asymmetrical channel if and only if A(x, y) > r 
for all X, y G C, X 7 ^ y. 

Our goal is to link these classical asymmetric codes to quantum AD codes. As discussed above, we start from 
the following definition. 

Definition 5: A classical code C over the alphabet B is called self-complementary if for any x G O, 10 x G C. 

For any self-complementary code C, there exists another code C such that C = {u+al: u G C, a = 0,... q— 1} 
and \C\ = q\C\. We may, for example, choose all u G C such that the first digit is 0. From C we derive the quantum 
code Q spanned by (I'i/'u): u S C} as given in Eq. (31). Our main result is given by the following theorem. 

Theorem 6: If C is a classical (linear or non-linear) self-complementary code correcting a single error with 
respect to Definition 3, then Q spanned by {|V’u): u G C} is a single-error-correcting code for the qudit AD 
channels 5 and A. 

Proof: Let Eij = \i){j\ with i,j G {0,1,..., g — 1} and i < j. For a small decay time r, in order to improve 
the fidelity of the transmission through the qudit AD channel A given by Eq. (4) from 1 — 0(r) to 1 — it is 

sufficient to correct a single £'i_i+i-error and detect one Z^-error for I = 1,2,... ,q — 1. The self-complementary 
form of \tjjf) given in Eq. (31) implies that A|'0u) = |'0u)- In turn, this implies that ('i/'vl^qlV'u) = 0 for any u, v 
and / = 1, 2,..., g — 1, i.e., the error-detection condition for a single Z* error is fulfilled. 

Next consider a single operator Every state of the quantum code is a linear combination of states |c) 

with c G C. Applying the operator Ei^i+i to |c) corresponds to a single asymmetric error. As the classical code 
C corrects a single asymmetric error, the distance A(u,v) between any two codewords u and v is at least two. 
Therefore, the supports (set of basis states with non-zero coefficient in the superposition) of the states l-ipf) and 
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£'j-*i+il'0v) are disjoint for all positions a, where denotes the operator acting at position a. Hence 

those states are mutually orthogonal. Finally note that for errors Ei^i+i acting on the same position, the operator 
ts diagonal and hence in the span of the operators which can be detected. I 

Corollary 7 : If there exists an {n,K, 3 )q self-complementary code C, then there exists an {{n,K/q))q quantum 
code correcting a single AD error. 

Such codes have, e.g., been studied in [14]. For linear codes, we have the following corollary. 

Corollary 8: If there exists an [n, k + 1, 3]g linear code C containing the all-one-vector 1 S C, then there exists 
an |n, fcjg CSS code correcting a single AD error. 

In the preceding corollaries we have used the notation (n, iF, d)q for a classical code of length n with K codewords 
and minimum distance d over an alphabet with q elements, and the notation [n, k, d]q = {n,K = q^, d) when the 
code is linear. 

VI. Single-Error-Correcting Codes: Even Lengths 

We now use Theorem 6 to construct some families of good single-error-correcting AD codes. Eor this, we 
need to find some good self-complementary single-error-correcting classical asymmetric codes. The best known 
direct construction of single-error-correcting codes for the binary asymmetric channel is the so-called Varshamov- 
Tenengolts (VT)-Constantin-Rao (CR) code [10, 39], with a natural generalization to g > 2. These VT-CR codes 
are non-linear codes, in both the binary and non-binary cases. 

Eor the binary case, many of these VT-CR codes are indeed self-complementary, and so they can be used to 
construct families of good single-error-correcting quantum AD codes [36]. As the VT-CR codes are nonlinear, the 
corresponding quantum codes are nonadditive codes. Unfortunately, for the non-binary case, the VT-CR codes are 
no longer self-complementary, so one needs some other constructions of good single-error-correcting for asymmetric 
channels. 

We will use the idea of generalized concatenation, which has been discussed in the context of constructing binary 
AD codes in [36], and in the context of constructing (classical) asymmetric codes in [18]. This method will allow 
us to construct good self-complementary asymmetric linear codes for the non-binary case, which will lead to good 
single-error-correcting quantum codes for AD channels. 

A. Qutrit Codes 

Eirst, we consider the case of g = 3. Eor the generalized concatenation construction, we choose the outer code 
as some ternary classical code over the alphabet {0,1, 2}, and the inner codes as: 

Co = {00,11,22}, Cl = (01,12,20}, Cj = (02,10,21}. (32) 

Then we have the following result. 

Theorem 9: Eor n even, generalized concatenation with an outer [n/2, fc, dja code results in an [n,n/2 -|- k]^ 
self-complementary linear code C. This code leads to an |n,n/2 + k— Ija quantum stabilizer code Q, correcting 
a single error for the channels S and A. 
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Proof: Note that Cq, Cj, and Cg are all self-complementary codes correcting a single asymmetric error. 
Therefore, any outer ternary code will lead to a self-complementary ternary code C, and hence a quantum code Q. 
A single amplitude damping error induces only a single error with respect to 0, 1, 2. As the outer ternary code has 
distance 3, such an error can be corrected. ■ 

Note that with respect to the symbols 0,1, 2, the induced channel is nothing but the ternary symmetric channel 


shown in Fig. 1. 



Fig. 1. The induced channel 7^3 for q = S (which is just the ternary symmetric channel), the induced channel 'IZ 4 for q = A, and the induced 
channel 7^5 for q = 6 . The arrows indicate the possible transitions between symbols. 


Example 10: For n = 6, take the outer code of length n/2 = 3 as {000,111,222} with distance 3. Generalized 
concatenation yields a self-complementary ternary linear code of dimension 4. The corresponding quantum code Q 
encodes 6/2 -f 1 — 1 = 3 qutrits. Both the best corresponding single-error-correcting quantum code |6, 2, 3 J 3 and 
the best possible asymmetric CSS code |6, 2, (3, 2}]3 (see Corollary 8) encode only 2 qutrits. 


B. The Case g > 3 

For g = 4, we choose the inner codes as 

Cg = (00,11,22,33}, Cl = (01,12,23,30}, 

Cj = {02,13,20,31}, Cg = {03,10,21,32}. (33) 

Similar as in Theorem 9, an outer code with distance three yields a self-complementary code from which a quantum 
AD code can be derived. Flowever, in this case, the induced channel for the outer code is no longer symmetric. A 
single damping error will, for example, never map a codeword of the inner code Cg to a codeword of Cg. So on 
the level of the outer code, there are no transitions between 0 and 2, or between 1 and 3. The induced quaternary 
channel 72.4 is shown in Fig. 1, where we see that errors only happen between ‘neighbors.’ 

The above constructions for g = 3,4 have a direct generalization to general g > 2. For a given g, choose the 
outer code as some code over the alphabet {6, 1,..., g — 1}. The g inner codes Cg, Cj,..., C-~i are the double¬ 
repetition code Cg = {00,11, ..., (g — l)(g — 1)} and all its g — 1 cosets Cj = Cg 0 (Oz), i.e., we apply the rule 
that Oz € Cp It is straightforward to check that each inner code has asymmetric distance 2, hence corrects a single 
asymmetric error. Similar as in the case of g = 4, a single damping error will only drive transitions between z,j 
for i = j ±\. For instance, for g = 5, the induced channel 72.5 is shown in Fig. 1. In general, we will write the 
induced channel as 72^ for outer codes over {0,i,...,g — 1}. 

Similar as Theorem 9, in general we have the following theorem. 
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Theorem 11: For n even, an outer [n/2,k]q code correcting a single error for the channel TZq leads to an 
[n,n/2 + k]q self-complementary linear code C and hence an |n,n/2 + k — 1], quantum code Q, correcting a 
single error for the qudit AD channels S and A. 

Note that the channel IZq is no longer a symmetric channel, so outer codes of Hamming distance 3 are no longer 
expected to give the best codes. It turns out, however, that single-error-correcting codes for the channel IZq are 
equivalent to single-symmetric-error correcting codes in Lee metric [5] (see also [24]), for which optimal linear 
codes are known (for a more detailed discussion, see [18]). 

VII. Single-Error-Correcting Codes: Odd Lengths 

The construction of AD codes for even lengths given in Sec. VI based on generalized concatenation is relatively 

straightforward. The inner codes are just 1-codes of length 2 with q codewords and their cosets. In [18], codes 

of odd length were obtained using a mixed-alphabet code, treating one position differently. This does not directly 
translate to the situation considered here, as the resulting code has to be self-complementary. 

Instead, we will use different inner codes, one of odd lengths and the length-two code from above. In particular, 
we can directly search for q mutually disjoint inner codes of length 3 which are 1-codes. 

For q = 4, consider the following Z 4 -linear code Co' of length 3 generated by {111, 002,020}: 

000 111 222 333 002 113 220 331 

020 131 202 313 022 133 200 311. (34) 

The code Cq' has asymmetric distance 2, as well as the three cosets Cy = Cq' -b 001, C 2 ' = Cq' + 010, and 
C'a' = C'o' +100. Applying generalized concatenation to the outer code (000', 111', 222', 333'} and the inner codes 
of length 2 and 3 for the first two and the third position, respectively, yields a self-complementary 1-code [7, 5 ) 4 . 
The corresponding quantum code has parameters |7,4]4. 

Note that the induced channel on the alphabet {O', 1', 2', 3'} is no longer 77 . 4 , but the symmetric channel over 
Z 4 . Therefore we have the following theorem for g = 4. 

Theorem 12: For n odd, an outer [{n — l)/2, fc, 3]4 code leads to an [n, (n + l)/2 + k]^ self-complementary 
linear 1-code C. The resulting quantum code Q = |n, (n — l)/2 -b fc ]4 corrects a single error for the qudit AD 
channels F and A. 

Proof: The inner codes Q of length two as well as the inner codes Cy of length three are self-complementary 
1-codes. The outer code has distance 3 which ensures that a single error mixing the inner codes can be corrected. 
For the outer code, we always take the last coordinate to be of type s', and all the other coordinates to be of type 
s, for s = 0,..., 3. Hence, the inner code for the last coordinate of the outer code has length three, while the other 
inner codes have length two. Therefore, for an outer [{n— l)/2, k, 3]4 linear code, generalized concatenation yields 
an [n, (n -b l )/2 -b fc ]4 self-complementary linear 1 -code C, corresponding to an |n, (n — l )/2 -b fcj 4 quantum code. 


We emphasize that the construction related to Theorem 12 is valid only for g = 4. For g > 5, the Zq-linear code 
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Co generated by {111, 013} and its q cosets are all self-complementary codes with asymmetric distance 2. For this, 
note that A(x,y) = 1 if and only if, up to permutation, x — y G {(1,0,0),(1,—1,0)}. For q > 5, the code Cq 
does not contain such a vector. Hence we obtain the analogue result as in Theorem 12 for q > 5. 

For <7 = 3 and q = 5, however, we cannot partition the trivial code [3,3]g into q self-complementary codes [3, 2]q 
with asymmetric distance 2, substantiated by exhaustive search. However, we can use an inner code of length five, 
resulting in the following theorem. 

Theorem 13: For q = 3,5 and n odd, an outer [(n — 3)/2, k]q code correcting a single error for the symmetric 
channel leads to an [n, (n -|- l)/2 -f k]g self-complementary code C and hence an |n, (n — l)/2 -f kjg quantum code 
Q, correcting a single error for the qudit AD channel A for q = 3,5. 

Proof: We can map the first digit of the outer code to q groups of codes of length 5. Again by exhaustive 
search, we find that for q = 3,5, we cannot partition [5, 5]^ into q self-complementary codes [5,4]^ with asymmetric 
distance 2, but we can get q codes [5,3],. For q = 3, the code Cq' is generated by (00011,01201,11111}, while 
for <7 = 5 it is generated by (00011,00102,11111}. 

Then our construction can be described as follows. For an [(n — 3)/2, k]q outer code, we use the length-5 code 
described above as the inner code for the first digit, and use length-2 code for the remaining (n — 5)/2 digits, 
leading to a code with length 1x5-1- x 2 = n. Similar as in the proof of Theorem 12, it follows that the 
resulting code corrects a single AD error. ■ 

In the nonlinear case, we can find larger codes. The results are summarized in Theorems 14, 15, and 16. 

Theorem 14: For q = 3 and n odd, an outer ((n — 3)/2,iT, 3)^ code leads to an (n,33 x self¬ 
complementary code. The resulting ((n, 11 x quantum code corrects one error for the AD channels 

S and A. 

Proof: An exhaustive search reveals that for g = 3, we can find three disjoint self-complementary codes of 
length 5 and asymmetric distance 2 with at most 33 codewords. Let 

Cq, = ( 00000 , 00011 , 00112 , 00220 , 01021 , 01110 , 

01202,02022, 02101,02120, 02211} (35) 

Then Cq' = (u -f al | u e Cosct = 0,1,2}, Cv = (u -f 00001 | u e Cq/}, C2' = (u -f 00002 | u e Cq/}. We 
construct the code similarly as in Theorem 13, i.e,. the inner code for the fist digit is C^/, and for the remaining 
(n — 5)/2 digits the inner code is C}. ■ 

For g = 5, we consider two constructions based on nonlinear codes. 

Theorem 15: Then for q = 5 and n odd, an outer ((n — l)/2,Ar, 3)^ code leads to an (n,20 X 5("-3)/2 a:)5 
self-complementary code. The resulting ((n, 4 x quantum code corrects one error for the AD channel 
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Proof: Although we cannot partition [3, Sjs, a weaker result can be obtained. Let 
Cq, = {000,002,020,022}, Cv = (001,004,021,024}, 

C2- = (003, on, 031,033}, Cz' = (010,023,041,043}, 

= {012,014,032,034}. 

Furthermore, set Ci' = {u + al | u e , a = 0,1, 2, 3,4}. 

Then we have a construction similar to Theorem 12, i.e., we use one copy of the length-3 inner code and (n—3)/2 
copies of the length-2 inner code. ■ 

Theorem 16: For q = 5 and n odd, an outer ((n — 3)/2,Ar, 3)5 code leads to an (n, 295 x self¬ 
complementary code. The resulting ((n, 59 x quantum code can correct one error for the AD channel 

n. 

Proof: In this construction we use self-complementary codes of length 5. By non-exhaustive search, we can 
hnd a self-complementary code with 295 codewords. The 59 codewords of Cq are shown in Table I. From those 
we derive Cq/ = {u -f al | u G Cq', a = 0,1, 2, 3, 4} and Ci' = {u -f OOOOi | u G Cq'}, 1 < i < 5. Then we have 

TABLE I 

Code Construction for 9 = 5 


00000 

00002 

00013 

00020 

00031 

00033 

00044 

00111 

00114 

00122 

00141 

00200 


00202 

00220 

00223 

00244 

00303 

00311 

00314 

00330 

00332 

00424 

00442 

01133 


01241 

01404 

01412 

02000 

02002 

02013 

02021 

02032 

02034 

02114 

02130 

02143 


02200 

02203 

02211 

02223 

02314 

02321 

02332 

02424 

02440 

03041 

03044 

03102 


03110 

03212 

03231 

03233 

03300 

03303 

03342 

03410 

03412 

03431 

04234 


a construction similar to that in Theorem 13. ■ 

As shown in Table II, for many lengths, the construction based on Theorems 9, 13, and 14 outperforms both the 
best known quantum codes with distance 3, and the CSS codes of Corollary 8. The dimension of the asymmetric 
quantum codes (AQECC) is taken from [14]. Only when n = 13, our construction performs worse than CSS and 
AQECC. This might be due to the fact that the outer code we use here is the CSS code [5,1, {3, 2}], which is not 
efficient since we also have CSS code |4,1, {3, 2}]. 

VIII. Multi-Error-Correcting Codes 

Eor the binary case, multi-error-correcting amplitude damping codes are discussed in [11]. The basic idea is that 
for the encoding |0 l) = |01), |1l) = |10), the amplitude damping channel simulates a binary erasure channel. 
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TABLE II 

DIMENSION OF SINGLE-ERROR-CORRECTING QUANTUM AD CODES FROM THE GF(3^) CONSTRUCTION WITH DISTANCE 3, THE CSS 
CONSTRUCTION, ASYMMETRIC QUANTUM CODES (AQECC), AND THE GENERALIZED CONCATENATION CONSTRUCTION (GC). 


n 

GF{3'^) 

CSS 

AQECC 

GC (linear) 

GC (nonlinear) 

4 

30 

3 ® 

1 

3 

3 

5 

3I 

3I 

6 

32 

11 

6 

32 

39 

11 

33 

33 

7 

33 

33 

29 

33 

11 X 3 ® 

8 

3 ^ 

3 '‘ 

84 

3® 

3 ® 

9 

35 

3 ® 

3 ® 

3® 

11 X 33 

10 

3 ® 

3 ® 

3 ® 

3® 

3 ® 

11 

3 ® 

3 ^ 

3 ^ 

3^ 

11 X 3 ® 

12 

37 

33 

38 

38 

38 

13 

38 

39 

39 

38 

11 X 3 ® 

14 

39 

39 

39 

310 

3I0 

15 

3I0 

3I0 

3I0 

310 

11 X 38 

16 

3II 

3II 

3II 

312 

312 


So one can use erasure-correcting code as outer codes to build codes correcting amplitude damping errors. In this 
section we consider generalizations of this construction, for both the binary and non-binary cases. 

It is mentioned in [II] that a possible generalization is to use |001), |010), |100) as the inner code, and a distance 
f + 1 quantum code as an outer code. However, it turns out that one can actually use |001), |010), |100), |111) as 
the inner code, and a distance f + 1 quantum code as an outer code. Here a single damping error will be treated as 
an erasure, and two damping errors or no damping can be treated as an error which is taken care of by the outer 
code. 

To generalize this idea to the case q > 2, one can take a similar approach. In this section we consider the 
channel A with Kraus operators given by Eq. (4). For g = 3, one can take the encoding |0 l) = |11), |1l) = |02), 
\2l) = |20). Then the amplitude damping channel simulates a ternary erasure channel. So one can use erasure- 
correcting codes as outer codes to build codes correcting amplitude damping errors. Actually, using a similar idea 
as the construction based on |001), |010), |100), |111) for the binary case, one can use |0 l) = |00), |1l) = |20), 
|2l) = |ll)j |3 l) = |02), and = |22) as the inner code, and still a distance t + 1 quantum code as an outer 
code. Here we give the general construction. 

Assume the length of the inner code is m. We choose the set 

S = {\a 1 a 2 ■ ■ ■ am) I fli + 02 + ... + Om is even} (36) 

as the orthonormal basis of the code, and let 

K=\S\. (37) 

If we have an outer code |n, k,t + l]if, we get an Inm, code correcting t errors for A. For g = 3, the code 
also corrects errors for the qutrit channel 5. We have the following theorem: 
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Theorem 17: The code constructed above corrects t errors for the channels A and S. 

Proof: We want to prove that the condition in Theorem 1 holds. Let Q be the outer |n, fc, f code. Take 

any two error operators and and for any two vectors li/ij and |'0j), we consider (^i|EjE|'|^j). Suppose 
that on some inner code, En has an odd number of errors while E^/ has an even number of errors. Then E|/|i/ij) 
will be in the space spanned by S while E^ji/ii) will be in the space perpendicular to S. So ('0i|E|E|/li/ij) = 0 
and the condition in Theorem 1 automatically holds. By symmetry, this argument also holds when Ei is even and 
Ef/ is odd. Now assume that on each inner code, the number of errors corresponding to E^ and E^/ have the same 
parity. We consider the series expansion of the operator E^Egr with respect to r. For this, we expand the tensor 
factor acting on the Eth particle such that Aij corresponds to the term of order rf Combined, we get 

E^El' = (Aiq + AiiT + yli2T^ + . . .) 

{A 20 + A 21 T + ^22'!'^ + • . •) 

0 ... 0 (AmO + AmlT + Am2T^ + . . .) (38) 

Note that each of Aiq, A 20 , ■ ■ ■ Amo is either the identity operator or the zero operator, which can be proved if we 
take the limit r —0. For the terms of the form Tt{A), there are at most t non-identity terms on the inner codes, 
which will become 0 according to the error-detection criterion. Thus the condition in Theorem 1 is satisfied. 

■ 

Finally we give an explicit expression for the dimension of the code K defined in (37) in terms of q and m. 

Theorem 18: 

\ /2 if q is even. 

K={ (39) 

[ (g'"-fl)/2 if g is odd. 

Proof: When q is even, we map \a 1 a 2 ... a^) to |g — 1 — ai, 02 ... Om), which is a one-to-one mapping from 
S to S. When q is odd, let i = min{j: aj f 0}, and we map |aia 2 ... a^f) to |ai... ai_i, q — a^, ... a^, 

which is a one-to-one mapping from S' \ {|00... 0)} to S. So |S| = |S| when q is even, and |S| = |S| -f 1 when q 
is odd. Thus |S| = g™/2 when q is even, and |S| = (g"* -f l)/2 when g is odd. ■ 

In Table III, we compare our codes with both stabilizer codes and AQECCs. We use the construction with g = 3 
and m = 2, and compare the codes that corrects t errors. From |n/2, fc, f -f IJs codes given (see also [29]), we can 
construct our ((n, K))^ codes where K = 5^', and we compare them with stabilizer codes |n, E, 2t + Ijja that have 
dimension K' = 3^ . It can be seen that our construction outperforms stabilizer codes, and our performance gets 
better in comparison for larger t. In most cases, our codes are better also better than the asymmetric CSS codes. 
Some of the CSS codes are optimal codes taken from [13], and others are based on the best known classical ternary 
codes [19]. For comparison, we also list an upper bound on the maximal dimension of an AQECC based on 
the known bounds for classical codes. In most cases, this bound cannot be achieved. 


TABLE III 

Dimension of quantum codes from our construction that correct t AD errors compared with stabilizer codes 
[n, k', 2t + IJs AND PURE ASYMMETRIC CSS CODES WITH PARAMETERS |n, k", {2t + l,t + 1}]3. 
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t 

n 

K 

loggX 

k' 

k” 

k" 

'^max 

2 

to 

5 

1.465 

1 

I10,1,{5,3}13 

2 

2 

12 

25 

2.930 

2 

112,3,{5,3}l3 

3 

2 

14 

125 

4.395 

4 

I14,4,{5,3}]3 

4 

2 

16 

625 

5.860 

5 

[16, 5, {5, 3)13 

5 

2 

18 

3125 

7.325 

6 

[18,7,{5,3}l3 

7 

2 

20 

15625 

8.790 

8 

[20,9,{5,3}|3 

9 

3 

14 

5 

1.465 

N/A 

I14,0,{7,4}]3 

0 

3 

16 

25 

2.930 

0 

[16,1,{7,4}|3 

1 

3 

18 

125 

4.395 

1 

I18,3,{7,4}]3 

3 

3 

20 

625 

5.860 

3 

[20,4,{7,4}|3 

5 

3 

22 

3125 

7.325 

4 

[22,6,{7,4}|3 

6 

3 

24 

15625 

8.790 

6 

[24,6,{7,4}l3 

8 

4 

18 

5 

1.465 

N/A 

N/A 

- 

4 

20 

25 

2.930 

N/A 

I20,0,{9,5}]3 

1 

4 

24 

625 

5.860 

1 

I24,4,{9,5}13 

5 

4 

26 

3125 

7.325 

2 

I26,4,{9,5}]3 

6 

4 

28 

15625 

8.790 

3 

I28,5,{9,5}l3 

8 

5 

26 

5 

1.465 

N/A 

I26,0,{13,6}13 

3 

5 

28 

25 

2.930 

0 

[28,1,{11,6}13 

5 

5 

30 

125 

4.395 

1 

I30,2,{11,6}13 

6 

5 

32 

625 

5.860 

1 

[32,4,{ll,6}l3 

7 

5 

38 

3125 

7.325 

4 

[38,7,{ll,6}l3 

12 

5 

40 

15625 

8.790 

6 

[40,8,{12,6}l3 

14 

6 

30 

5 

1.465 

N/A 

7 

0 

6 

36 

25 

2.930 

0 

7 

6 

6 

38 

125 

4.395 

1 

I38,2,{13,7}13 

8 

6 

40 

625 

5.860 

1 

I40,3,{13,7}13 

10 

6 

42 

3125 

7.325 

2 

[42,5,{13,7}l3 

11 

6 

44 

15625 

8.790 

2 

[44,6,{13,7}l3 

13 


IX. Applications to other AD channels 

Our method can also be applied to some other amplitude damping processes. For instance, when g = 3, in 
addition to the channel S, there are the following natural decay processes of three-level atoms; the A-pattern or the 
V pattern. They are illustrated together with the S channel in Fig.2; 

For the A-pattern, the master equation is 

— = ki(2a^2P^12 ~ '^12'^12P ~ P^12^12) 

+ k2{2aQ2P<^02 ~ O'02^02P ~ P^02'^02)^ 









| 2 > | 2 > 
| 1 > 

| 0 > 
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-| 1 > 


| 2 > 


-| 0 ) 


| 1 > 


S V 

Fig. 2. Decay processes for three-level atoms with different level-structures. 


-| 0 > 


where ai 2 and a ^2 are the same as above, and 

ao -2 = |0)(2|, a+ = |2)(0|. (40) 

By direct calculation, one can verify that the evolution of this master equation can be expressed by using the 
following Kraus operators: 

2 

P('^) = (41) 


2=0 


where 


^0 = diag{l, 1 , \/i-Ti-^}> 

Ai = v^| 0 )( 2 |, 

A 2 = v^|l)( 2 |, 


and 


71 = 

72 = 


^2 i_ g-2(fci+fe2)r 
I _ g-2(fcl+fe2)T 


fci + fc 2 - 

ki 


= 2k2T + 0{t^), 
= 2kiT + 0 {t^). 


ki + k2 - 

Both 7 i and 72 are of first order in r. 

For the V-pattern, the Kraus expression has been found in [ 8 ], which is given as 




2=0 


where 


^0 = diag{l, \/r^^}, 

Ai = yyr|o)(i|, 

A 2 = V^|0)(2|, 


and 


71 = 1 - = 2kiT + O(t^), 

72 = 1 - = 2k2T + O(t^). 

Again, both 71 and 72 are of first order r. 


(42) 


(43) 


We also introduce the classical counterparts of these two channels which are shown in Fig. 3. Here the arrows 
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indicate the allowed transitions. The classical channels Li and £2 correspond to the amplitude damping channels 
V and A, respectively. 


1 



2 



0 


C2 


Fig. 3. Classical channels for trits: the arrows indicate the allowed transitions. 


For the l/-channel and the A-channel, we construct codes using similar ideas as shown above. For single-error- 
correcting codes, we can still use the idea of self-complementary codes based on the corresponding classical codes, 
which have been studied in [33] and [22]. Unfortunately, none of these constructions can be adapted to make 
self-complementary codes. For codes with short lengths, one can use numerical search. 

For multi-error-correcting codes, we want to choose the basis of the inner code so that one single damping error 
will project the state to an orthogonal subspace. Let Tq and Ti be the set of binary strings of length m with even 
and odd parity, respectively: 

Tj = {xiX2 ... e {0,1}™ \ xi ® x^ ® ® Xm = i}- (44) 

For the channel £ 1 , a binary 0 is mapped to the channel symbols 1 and 2, and a binary 1 is mapped to 0. For the 
channel £ 2 , a binary 0 is mapped to the channel symbols 0 and 1, and a binary 1 is mapped to 2. Then Tq and 
Ti are mapped to the set of codes Sq and which are two sets of possible inner codes. 

Example 19: For m = 2, we have 

To = {00,11}, Ti = {01,10}. (45) 

For the channel £ 1 , we get 

So = {00,11,12,21,22}, Si = {01,02,10,20}, (46) 

and for the channel £2 channel: 

So = {00,01,10,11,22}, Si= {01,02,10,20} (47) 

We have the following theorem. 

Theorem 20: For any codeword w S 5^, i = 0,1, when an error of the corresponding channel occurs, the 
resulting string v will be in Si-i. 

Proof: We only prove the case when the channel is £ 1 , the proof for the channel £2 follows using the same 
argument. Let w = W 1 W 2 ■ ■ ■ Wm and v = V 1 V 2 ■ ■ - Vm, and let the error occur at position t. Then wt is either 1 
or 2, while vt = 0, and Wk = Vk, for all k t. Suppose that in our construction, wi ... Wm corresponds to the 
binary string oi... Om, and vi .. .Vm corresponds to the binary string bi... b^, then we have at = 0, bt = 1, and 
CLk = bk, for k Since a G Ti, we know that b G Ti_i, and hence v G S'l-i. ■ 

Corollary 21: For any codeword w G 5^, i = 0,1, when an odd number of errors of the corresponding channel 
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occur, the resulting string v will be in Si-i; when an even number of errors of the corresponding channel occur, 
the resulting string v will be in Si. 

With arguments similar as those in the proof of Theorem 17, for i = 0,1, if we use the quantum code Qi spanned 
by {|u): u € Si} as the inner code and an |n, k,t + ll|Si| code as the outer code, we get an |nm, code that 
corrects t errors for the corresponding quantum AD channel. 

Finally we will give an explicit expression for |S'i| in terms of m. To show its dependency on m, we write 
CXm — 1*^01 Ultd 
Theorem 22: 

am = |5'o| = i(3'" + l) and = |,5i| = i(3”^ - 1). (48) 

Proof: For any string a = 0102 ... Om & To, either am = 0 or am = 1- For am = 0, the strings in So that a 

maps to are just the strings oi ... Om-i maps to, concatenated with a single symbol chosen from two options {1^2 

for Cl, and 0,1 for £ 2 )- For am = 1, the strings a maps to are those that ai ... am-i maps to, concatenated with 
a fixed symbol. So we have the recurrence relation 

— 2.CXm—l T /Tm —1- (49) 

We also have 

am+/3m = 3™. (50) 

Solving Eqs. (49) and (50) together with the initial condition ai = 1 and /3i = 2, we have 

am = ^(3™ + !) and = 1(3™ - 1), (51) 

which proves the theorem. ■ 

To achieve the maximal size of the constructed code, we should always choose {|u): u € S'o} as the inner code. 
Example 23: We take m = 2, and the outer code is [5,1, Sjs (see, e. g., [7]), which is 

1 . ^ 

|fc) —-= + 9 + fc) (g) Ip + r) (g) |g + r) 0 Ip) (g) |g), (52) 

p,9,r=0 

where w = exp(27ri/5). We substitute |0), |1),..., |4) with |00), |11), |12), |21), and |22), respectively, and get a 
((10, 5))3 code which corrects 2 errors of the channel V. If we substitute |0), |1),..., |4) with |00), |01), |10), |11), 
and |22), we get a ((10, 5))3 2-code for the channel A. In comparison, the best stabilizer code of length 10 which 
corrects 2 errors is |10,l,5j3. 

Other examples for which our construction outperforms stabilizer codes include the ((12, 25))3 and ((14,125))3 
2-codes constructed from outer codes |6,2,3j5 and |7, 3,3]5, while the best corresponding stabilizer codes are 
[12,2,5]3 and [14,4,513. 
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